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Abstract 

The Direct Simulation Monte Carl<- method is used to study the 
hypersonic, rarefied flow interference effects on a flat plate caused by 
nearby surfaces. Calculations focus >>n shock-boundary-layer and shock-lip 
interactions in hypersonic inlets. Results art presented for geometries 
consisting of a flat plate with different leading-edge shapes over a flat 
lower wall and a blunt-edged flat plate over a 5-degree wedge. The problems 
simulated correspond to a typical entry flight condition of 7.5 km/s at 
altitudes of 75 to 90 km. The results show increases in predicted local 
heating rates for shock-boundary-layer and shork-lip interactions that are 
qualitatively similar to those observed experimentally at much higher 
densities. 


Nomenclature 

Cfl heat-transfer coefficient, 2q/p 00 U 00 

Cp pressure coefficient, 2p/ ooU«,2 

Knco freestream Knudsen number, A«>/Rn 

£ characteristic dimension 

Moo freestream Mach number 

p pressure 

q heat flux 

Rn leading-edge nose radius 

Re,*, Reynolds number, pooUoo^/ poo 


T 


temperature 


Uoo freestream velocity 

Xi mole fraction of species i 

x coordinate measured from leading edge of flat plate 

Xoo freestream mean free path 

Poo freestream viscosity 

Poo freestream density 

0 angular position on circular leading edge 

Subscripts 

o stagnation value 

is isolated value 


Introduction 

Viscous-inviscid interactions occurring in high Mach number flight 
represent one of the most severe problems associated with hypersonic cruise 
and entry vehicles. The strongest interactions are usually associated with 
some form of shock impingement such as that occurring on the cowl lip or wall 
boundary layers of hypersonic inlets as illustrated in Fig, 1, These 
interactions can result in high local heating which may lead to component 
failures. An excellent review of the variety of problems that can occur 
along with examples of some actual failures is given by Korkegi, 1 

The most severe conditions for vi scous-inviscid interaction problems are 
usually associated with the higher densities of the continuum flow regime. 
However, similar interactions are expected to occur at the lower densities 
associated with initial entry or sub-c rbital maneuver conditions (typically 
above 70 km). Although the strength of the interactions may be expected to 




be weaker in these transitional flows, localize! problems may still occur. 
Furthermore, even when the overall flovfield is essentially continuum, the 
analysis of the problem may require modifications to account for local 
transitional or rarefied flow effects. » The ase of a purely continuum- 
based analysis typically leads to an overprediction of the local heating 
rates in such cases. 

The Direct Simulation Monte Carlo (DSMC) m;thod of Bird 4 * 5 has been 
used extensively to study a variety of transitional and rarefied flow 
problems in hypersonic flows. For the most part, these studies have focused 
on modeling the physics 4-9 or on calculations of complex flows about simple 
geometries. 10-1 ^ Studies of more complex geometries have been limited since 
they require the more difficult three-dimensional codes or require prohibi- 
tively large computer storage and CPU time to resolve the flow around complex 
shapes even for two-dimensional problems. However, with the current availa- 
bility of supercomputers and with recent improvements to existing DSMC 
computer codes * that provide easier setup for complex geometries, 
multi-surface interaction problems can be more readily studied. 

The purpose of the present paper is to present initial results of a 
study of some of the basic interactions shown in Fig. 1 for hypersonic 
inlets. Although calculations are presented for specific geometries (see 
Fig. 2), the results are mainly intended to demonstrate the nature of the 
interactions at low densities and to show both the similarities and the 
differences in comparisons to analogous interactions in continuum flows. A 
flat plate was used to simulate an inlet cowl, since extensive studies have 

a Bird , G.^A.T" Private Communication, !987. 
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been made of the flow over isolated flat plates and the expected flow 

12 13 1 fe 21 

behavior is reasonably well understood. > > 

The results in this paper focus on some of the basic interference 
effects between two parallel flat plates and on the interference caused by an 
oblique shock impinging on the leading edge of the plate. The calculations 
were performed at nominal entry-type conditions of 7.5 km/s at an altitude of 
90 km for the parallel flat plates and over an altitude range of 75 to 90 km 
for the shock impingement calculations. 

Direct Simulation Monte Carlo Method 
The DSMC method of Bird has been described extensively in the 
literature. 4 * 5 The method consists of simulating a gas by some thousands 
of molecules whose positions, velocities, and internal states are stored and 
modified with time as the molecules move, collide, and undergo boundary 
interactions in simulated physical space. A time step is chosen such that 
the movement and collision processes can be assumed to be uncoupled over the 
duration of the step. The calculations are started with an initial state 
such as a vacuum or uniform equilibrium flow. The simulation time is 
identified with physical time in the real flow, and all calculations are 
treated as unsteady. When boundary conditions appropriate to steady flow are 
imposed, the solution is the asymptotic limit of the unsteady flow. 

A physical space computational cell network is required only to 
facilitate the choice of collision pairs and the sampling of flow proper- 
ties. The cells may be divided into rubcells for selection of collision 
pairs to better simulate the nearest neighbor collisions in regions of high 
vorticity. The collisions are simulated using the "variable hard sphere" 


(VHS) molecular model with modifications to account for rotational and 
vibrational states in diatomic or poly itomic gases. 5 

Geometries and Conditions for the Calculations 

Two basic two-dimensional geometries, shown in Fig. 2, were studied. 
Figure 2(a) shows a sharp-edged flat plate parallel to a flat lower wall. 

This geometry was used to study the interactions between the lower wall and 
the flat plate with different boundary conditions imposed on the lower wall. 

The effects of the leading-edge shapes shown in Fig. 2(c) were also 

computed. A set of calculations was also made with the lower wall 
representing a 5-degree wedge ramp (Fig. 2(b)) to produce an oblique shock 
that impinged on the leading edge of the plate. For these calculations, only 
the blunt (circular) leading edge was used. 

Calculations for the flat lower wall were performed only for the 90-km 
condition using the freestream values given in Table 1. Calculations for the 
5-degree wedge lower wall were made at all four altitudes given in Table 1. 

The freestream Mach number was 27.2 in all casei. The species mole fractions 

at altitudes below 90 km were held constant at values appropriate to 90 km 
using the tables of Ref. 22. Thus the chemical species simulation is 
realistic only at the 90-km condition. Some of the flow parameters, 
expressed in terras of the leading-edge nose radius, are given in Table 2. 

The computational domains and typLcal cell grids for the two geometries 
are shown in Fig. 3. The heavy solid lines represent region boundaries 
within which the time step and the molecule weighting factor (relating the 
number of computational molecules to the number of physical molecules) are 
constant. For the circular leading-edge cases, cells were clustered around 
the leading edge. The cell size was reduced at the lower altitudes in order 


5 



to minimize the flow gradients across the cell in a direction normal to the 
leading-edge surface. 

For all calculations, the wall temperature of the flat plate was assumed 
to be constant at 1000 K, and gas-surface collisions were assumed to be 
diffuse with full thermal accommodation. The same wall conditions were 
generally used for the lower wall except for a few cases with the flat lower 
wall in which the boundary was treated as specular or as a freestream 
boundary. 


Results and Discussion 
Isolated Flat Plate Predictions 

To establish the undisturbed nature of the isolated flat plate flow, 
calculations were made at the 90-km conditions with each of the leading-edge 
shapes and with the lower wall treated as a freestream boundary. The 
pressure and heat transfer distributions along the upper surface of the plate 
for these cases are shown in Fig. 4. At 90 km, the length of the plate 
corresponds to about 25 Aoo, and the results are typical of transitional 
flow behavior. ’ For the -20 degree wedge, the upper surface is flat, 
and the results are nearly identical to the zero thickness case. The slight 
differences between the -20 degree and zero-thickness results are attributed 
to the fact that molecules striking the lower surface of the leading edge 
have a fairly high probability of reaching the upper surface through gas-gas 
collisions at this Knudsen number. With the +20 degree wedge and circular 
leading edges, the pressure and heat-transfer coefficients are highest near 
the leading edge and decrease downstream, approaching the zero thickness 


^Moss, J. N. , Unpublished Results, 1907 



results near the end of the plate. While the +20 degree wedge shows higher 
pressures and heat transfer in Fig. 4 than the circular leading edge, results 
are shown only for the flat portion of the plate downstream of the leading 
edge. The peak stagnation values on the leading edge itself are actually 
higher for the circular leading edge. In all cases, the pressure and heat- 
transfer coefficients are much higher than the calculated free-molecular 
values of 0.003 for Cp and 0.01 for Ch at these conditions. The trends 
shown in Fig. 4 are qualitatively similar to those expected in continuum 
flows. 

Interactions with Flat Lower Wall 

The pressure and heat-transfer distributions along the lower surface of 
the plate with the blunt leading edge are shown in Fig. 5 for different lower 
wall boundary conditions. The freestream lower boundary condition represents 
the isolated or undisturbed flat plate. A specular lower boundary is 
equivalent to a plane of symmetry, and the geometry for this calculation is 
analogous to the entrance of a two-dimensional channel whose height is twice 
the distance between the plate and the symmetry plane. The results show a 
gradual rise in the pressure and heat transfer due to the growth of the 
boundary layer along the lower surface of the (upper) flat plate. 

With a diffuse lower wall, the geometry becomes similar to an inlet with 
a flat ramp, and the results show a more complex behavior. The pressures are 
much higher than the undisturbed predictions, and the distribution shows a 
steep rise starting at about 0.16 m (approximately 7 Aco). The pressure 
reaches a maximum at about 0.32 m and then decreases continuously to the exit 
of the channel. The heat-transfer distribution exhibits a maximum near the 
leading edge and decreases continuously downstream. However, there does 
appear to be an inflection in the heat -transfer distribution near the 
location of the pressure rise. 
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The above behavior is caused by Che interaction of a diffuse bow shock 
generated by the leading edge which strikes and is reflected , from the lower 
wall. This reflection is evident in the pressure distribution along the 
specular lower wall shown in Fig. 6. Significant differences in the location 
of the reflected bow shock are seen for different leading-edge shapes. These 
differences are probably more pronounced at the 90-kra conditions than might 
be expected at lower altitudes because of the large curvature and thickness 
of the shock, which is roughly the same order of magnitude as the channel 
height. 

For the diffuse lower wall (see Fig. 7), the differences in the 
Interaction for different leading-edge shapes are not as pronounced. Also, 
the relative strength of the interaction (in terms of the ratio of pressures 
across the interaction region) is reduced by the presence of the lower wall 
boundary layer (although the overall pressures are higher due to 
boundary-layer displacement effects). The results shown here are qualita- 
tively similar to the experimental data of Kaufman and Johnson for weak 
shocks incident on laminar boundary layers. Although the experiments were 
conducted at much higher densities (Re^lO^ per meter) and at a lower 
Mach number (M»“8), a reasonable comparison can be made based on the peak 
pressure rise. For the results in Fig. 7, the ratio of peak to upstream 
pressures is about 8 to 10, while the ratio of the peak to upstream heating 
rates is about 2.5 to 3.0. Based on the correlation of heating-rate amplifi- 
cation to peak-pressure amplification given in Ref. 23, the heating-rate 
amplification at the same pressure amplification would be roughly 15 to 20 at 
the Reynolds numbers of the experiment. The weaker interaction for the 
present calculations is perhaps not too surprising, since the shock is much 
more diffuse at the lower densities. 
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The pressure and heat-transfer distributions downstream of the shock 
interaction on the diffuse lower wall shows significant streamwise gradients 
extending to the exit of the internal region. Since the boundary layers are 
quite thick relative to the channel height, a significant portion of the 
internal flow is subsonic. Therefore, it is quite likely that the location 
of the exit and the boundary conditions imposed on the exit plane may have a 
significant influence on the shock -boundary layer interaction. Because of 
this possible influence, the post-shock flow quantities may not be 
representative of those for isolated shock reflections from laminar boundary 
layers. 

Interaction with 5-Degree Wedge Lower Wall 

For the calculations with the 5-degree wedge lower wall, the main 
interest in this study is on the interaction of the oblique shock generated 
by the wedge with the leading edge of the flat plate. Therefore, the 
following results focus on the surface properties on the blunt lip of the 
plate. 

The pressure and heat-transfer distributions around the leading edge of 
an isolated flat plate are shown in Fig. 8 for several altitudes. The 
stagnation pressure coefficient shows a slight decrease with decreasing 
altitude, whereas the stagnation heat-transfer coefficient decreases somewhat 
more rapidly. (It should be noted that although Cjj decreases with 
decreasing altitude, the absolute heating rate actually increases.) The 
results shown here are in reasonable agreement with the predictions of Cuda 
and Moss 10 that were obtained under similar freestream conditions but with a 
larger nose radius. A correlation of the stagnation heat-transfer 
predictions for the present calculations with those of Cuda and Moss is 
presented later in this paper. 
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As a result of the viscous shock layer caused by the 5-degree wedge, the 
leading edge of the flat plate is subjected to higher temperatures and 
densities. The qualitative features of the flowfield are demonstrated by the 
gray-shaded contours shown in Fig. 9, where the darker regions represent 
higher temperatures and densities. The variation in the extent of the 
overall thermal environment with altitude is clearly shown by the temperature 
contours, while the thickness and location of the oblique wedge shock are 
reasonably evident in the density contours. At 90 km, the distance from the 
start of the 5-degree wedge to the leading edge of the flat plate corresponds 
to about 50 A co. At 75 km, this distance corresponds to about 800 Aeo. At 
90 km, the shock is very thick relative to the leading-edge radius, and the 
peak density in the undisturbed shock is about twice the freestream value. 

At 75 km, the shock is more clearly defined with a peak density that is about 
four times the freestream density. 

It should be noted that the gradients in flow quantities across the 
shock are probably not adequately resolved at the lower altitudes because of 
the relatively coarse grid used for sampling. Although an attempt was made 
to accurately simulate the collisions in the shock by increasing the number 
of subcells at the lower altitudes, the average flow properties are always 
determined from the total sample within a cell. Since the gray contours are 
produced by filling each cell with a level-of-gray halftone shading that is 
proportional to the flow quantity, the shock may appear to be artificially 
smeared in these contours. However, the cell size used near the leading-edge 
surface was adequate to give reasonably accurate predictions of surface 
heating as determined by a cell-size sensitivity study. 

The pressure and heating-rate distributions around the plate leading 
edge with the 5-degree wedge shock are shown in Fig. 10. Both quantities 
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have been nondimensionalized by the peak stagnation values for the isolated 
leading edge at the corresponding altitude. Therefore, peak values greater 
than one indicate an amplification of the pressure or heating rate resulting 
from impingement of the shock. At 90 and 85 km, there is little effect of 
the shock, while the 80- and 75-km results show significant amplification of 
both pressure and heating rate in the stagnation region. 

A comparison of the heat-transfer coefficients with the 5-degree wedge 
with those for the isolated leading edge is given in Fig. 11 as a function of 
the product of the freestream density and the leading-edge radius. The solid 
line is taken from the results of Cuda and Moss, 10 and the present results 
correlate well for isolated blunt leading edges. The impingement of the 
wedge shock causes the heat-transfer coefficient to increase significantly 
with increasing density. Since the actual heating rate is proportional to 
the product of Cjj and the freestream density, the increase in heating is 
quite significant. 

Shock impingement effects on circular leading edges have been studied 
extensively for hypersonic flows at much higher densities than the present 
study. The works of Edney, 21 * Johnson and Kaufman, 25 and more recently 
Wieting and Holden describe in detail the nature of the interaction at 
continuum flow conditions. Although a direct comparison of the flow details 
from the present DSMC calculations with those given in these references is 
inappropriate, a qualitative comparison of the correlation of heating-rate 
amplification with pressure amplification may be appropriate. Using such a 

correlation, a comparison of the DSMC results with the shock-boundary-layer 
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Interaction data of Kaufman and Johnson and the shock-lip Interaction data 
of Wieting and Holden is given in Fig. 12. Although no direct conclusions 
should be drawn from this comparison, the interaction effects predicted by 
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the DSMC method do show similar trends to those for shock interactions with 


laminar boundary layers. It also appears that the DSMC results adequately 
represent the shock-lip interactions !:hat might be expected at the lower 
densities. Further DSMC calculations at the higher densitites and/or 
experimental measurements at the lower densities are needed. 

Concluding Remarks 

This initial DSMC study of interference effects on the hypersonic, 
rarefied flow about flat plates has focused on viscous-inviscid interactions 
relevant to inlet-type geometries. The results have been interpreted from 
somewhat of a continuum viewpoint. However, for the flight conditions 
presented, the distinction between viscous and inviscid regions is difficult 
at best. Therefore, the DSMC results need to be examined in more detail to 
better understand the physics involved in "viscous-inviscid” interactions at 
low densities. 

The geometries studied and the boundary conditions used do not 
necessarily represent realistic hypersonic vehicle problems. This is 
particularly true for the internal flow simulations where the mass flow rate 
depends strongly on the internal geometry and exit boundary conditions. The 
physical scale of the geometry is also important, and calculations for much 
higher Reynolds numbers would be useful. However, the DSMC method is not 
currently practical for very large vehicles at the high Reynolds numbers. 
Therefore, methods which can combine the computational capabilities of the 
DSMC method with more conventional continuum methods are needed. 

In spite of the limitations mentioned, the results of this initial study 
have demonstrated the usefulness of the DSMC method for the analysis of 
complex interaction problems such as those present in hypersonic inlets. The 
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results show qualitative trends that appear to be similar to those given by 
available experimental data. 
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Table 1. Freestream conditions 



Altitude , 
km 

Density, 
kg/m ^ 

U , 

00 ’ 

kra/s 

T , 

oo' 

K 

Mole Fractions* 

Reji 

m - 1 

\ , 

oo' 

m 

X °2 

x N 2 

*o 

90 

3.4xl0 -6 

7.5 

188 

0.209 

0.788 


2.0x10 3 

2.4x10- 

85 

8.2xl0 -6 

7.5 

188 

0.209 

0.788 


4.9x10 3 

9.9x10" 

80 

1.8xl0" 5 

7.5 

188 

0.209 

0.788 

0.003 

l.OxlO 4 

4.4x10“ 

75 

4.0xl0“ 5 

7.5 

188 

0.209 

0.788 


4.3x10** 

2.2x10- 


2 

3 

3 

3 


*Held constant at values appropriate for 90 km. 


Table 2. Parameters based on nose 

radius of blunt-edged flat 
plate, 

[Rfl * 0,005 m] 


Altitude, 

km 

Re 

00 

Kn 

OO 

90 

10 

4.74 

85 

24 

1.98 

80 

52 

0.88 

75 

109 

0.41 
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Fig, l Typical shock interactions in hypersonic inlets. 
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Fig. 2 Concluded. 
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Fig. 3 Typical grids for DSMC calculations. 
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(b) Heat transfer distributions. 


Pig. 5 Concluded. 




Fig, 6 Interact Lon of leading-edge shock with a specular lower wall 
boundary. (Altitude - 90 km) 
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(a) Pressure distributions. 


Fig. 7 Interaction of leading-edge shock with a diffuse lower wall 
boundary. (Altitude = 90 km) 
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(b) Heat transfer distributions. 


Fig. 7 Concluded. 
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(a) Pressure distributions. 


Fig. 8 Isolated blunt leading edge. 
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(a) Temperature. 


Fig. 9 Gray-shaded contours of flowfield for 5-degree wedge lower wall. 
(Darker shades denote higher values.) 
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Fig. 9 Concluded. 
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Fig. 10 Interference effects of 5-degree wedge on blunt leading edge. 
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altitudes of 75 to 90 km. The results show increases in predicted local 
heating rates for shock -boundary-layer and shock-lip interactions that are 
qualitatively similar to those observed experimentally at much higher 
densities. 
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